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The charm quark's mass is determined from Monte Carlo calculations of the cc spectrum. The main sources 
of uncertainty are perturbation theory (for conversion to MS), the continuum- limit extrapolation, Monte Carlo 
statistics, and the effects of quenching. The (preliminary) result for the MS mass is m c h(m c h) = 1.33 ±0.08 GeV. 
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1. INTRODUCTION 

Nonperturbative lattice calculations of the ha- 
dron spectrum provide a connection between ex- 
perimentally measured masses and the couplings 
of the (lattice) QCD Lagrangian. By convention, 
however, the MS couplings a(Mz) and m(//), 
used in phenomenology, are usually quoted. The 
two sets of definitions can be related to each other 
in perturbation theory. For example, 



m(p) = M(a) 



1 



(C + 70 ln/i a ) 



4tt 



(1) 



where 70 = 4. The lattice mass M and, by impli- 
cation, Co = Cq(M) are specified below. Eq. ([!]) 
omits higher orders in the gauge coupling a and 
power-law artifacts. 

This paper determines the charm quark's mass, 
m c h, from quenched calculations of the cc spec- 
trum. To anticipate the main sources of uncer- 
tainty, let us recall recent determinations the av- 
erage of the up and down quarks' masses Jlj||. 
There the three largest uncertainties JjJ stem 
from, in descending order, the quenched approxi- 
mation, the extrapolation to the continuum limit 
(even with the clover action), and perturbation 
theory. Each of these takes on a different guise 
for charm, however. 

The error in a coupling from quenching can 
be partly explained by noting that couplings run 
differently in the quenched approximation 
One can account for this effect by running the 
couplings down to typical mesonic momenta with 
rif = and then back up to a high scale with 
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/ 0. But m(/i) does not run for /1 < m, so 
quenching should not affect m(m) much 

One might expect lattice spacing errors to be 
worse for charmonium than for light mesons, since 
0.4 < mo,ch a < f on our lattices. Our spec- 
trum calculations are of mass splittings and the 
so-called kinetic mass of the meson, for which 
the cutoff effects are powers of \pa\, not ra^a Q. 
Indeed, we exploit two methods for determining 
?Ti c h, with opposite cutoff dependence. The two 
continuum limits agree, so cutoff effects are under 
better control. 

That leaves perturbation theory as the source 
of the largest uncertainty. To make the most of 
the one-loop approximation, the only order avail- 
able, we use results for mod ^ @,|. Further- 
more, we try to reduce the effect of truncating at 
one loop by choosing a(q*) in Eq. (Q) to absorb 
logarithms from higher orders ||[l0) . 

2. CUTOFF EFFECTS 

In a heavy-quark system, such as charmo- 
nium, typical three-momenta are only a few hun- 
dred MeV, suggesting that worrisome lattice ar- 
tifacts are of order (m^a) s . On the other hand, 
it is well-known that actions for Wilson fermions 
approach the static limit as mod — > 00, showing 
that higher-dimension operators are suppressed 
by a factor of order l/(moa) r . The lattice Hamil- 
tonian (defined by the transfer matrix) clarifies 
the middle ground, mo a « 1. One finds [§ 

H\ at = H cont + SH. (2) 

Contributions to the artifact 5H take the form 

(a6Hm)~g 2l bm( mo a)\pa\ s " + \ (3) 
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where p is a few hundred MeV, and s n > 0. The 
function 6„ is bounded || • It is safe to replace it 
by a number of order unity, and thus the effect is 
about the same size for splittings of charmonium 
as for masses of light-quark hadrons. 

Eq. (H) applies only if the hopping parameter k 
is adjusted until the meson's kinetic mass 



from the meson's kinetic mass: 



(4) 



M 2 := (d 2 E/dp^l 

equals the meson's physical mass. When Ma ^ 0, 
the rest mass Mi :— E(0) is smaller. Never- 
theless, the splittings of meson rest masses are 
accurate up to Eq. (^). In particular, the spin- 
averaged binding energy 

) MC -2(M 1Q a) PT , 



B ia := (M 



(5) 



where M^q is the spin average of mesons' rest 
masses, has relative errors of order min(p 2 a 2 , v 2 ) 
p"l| ]. (When the quark's rest mass M\q is com- 
puted to finite order in perturbation theory, B±a 
suffers perturbative errors as well.) 

To determine m c h we rely, therefore, on the fol- 
lowing Monte Carlo calculations: We define the 
lattice spacing a in physical units from AM = 
M hr - \{M Vc + 3M J/4 ,) §. We then obtain the 
quark mass either from the spin-averaged binding 
energy B\ of the IS states, or from their spin- 
averaged kinetic mass M 2 qq. 

3. PERTURBATION THEORY 

3.1. When m a ^ 

If the Monte Carlo has mofl ^ it is necessary 
to take mod ^ when deriving Eq. ([!]). Although 
Co remains bounded its value can change 

significantly for nonzero Too a. 

Eq. (]]]) is obtained by computing the quark's 
pole mass in lattice and in MS perturbation the- 
ory. Because the lattice breaks Euclidean invari- 
ance, several "masses" (Mi, M 2 , etc) describe the 
pole. One would like to pick a pole mass without 
dire lattice artifacts. We use two methods. In the 
first, we take the binding energy and set 



TOp } c 



(G) 



QQ 

with B\a from Eq. (|J) and a from AM. In the 
second method, we use the quark's kinetic mass, 
but reduce uncertainty in tuning k by taking a 



m-polc = {M 2 qo) 



pt 



,^expt 

QQ 



( M 2QQ a ) 



MC • 



(7) 



When B\a and M-iga are expanded in perturba- 
tion theory, Eqs. (J6|) and (f7[) can be matched to 
the expansion of m po i e in MS. The manipulations 
at one loop define M and Co in Eq. (Q). 

One needs, therefore, the loop corrections to 
the quark's rest and kinetic masses. From formu- 
las j^H f° r Mi and M 2 , to all orders in g 2 and 
in toq a, one can expand 



M x => afMf. 



One finds m[ 0] = log(l+M ), where M = 1/2k- 
1/2/w- Refs. @,| show results for m[ 1] . 

The kinetic mass has further loop corrections, 
so let Z M2 (Mi) = M 2 /to 2 (Mi). The function 
to 2 (M) is chosen so that — 1 , but it is eval- 
uated at the all-orders Mi. With this definition 
Zm 2 {0) — 1, to all orders in g 2 . Also, Z^ 2 is 

tadpole-free. It is small (0 > Z^} 2 > —0.1) and 
hardly depends on the clover coupling csw f§- 

3.2. Choosing a(q*) 

With only the one-loop approximation at hand, 
the right-hand side of Eq. (0) is sensitive to the 
choice of scheme for a and its scale q* . Since 
Eq. (|l|) is the combination of lattice and MS per- 
turbation theory, the original series must be ex- 
pressed in a common schem and the scales must 
be run to a common one. Here we use the scales 
suggested in Refs. primarily for ay, but 

also for a. 

For dimensional regulators Ref. B] prescribes 



E 



(8) 



lni ?BLM/^ 2 = I* /I, 



(9) 



where /* is derived from the Feynman diagram 
for I by replacing gluon propagators by 



\n(q 2 /p 2 ) - b{/pl 



(10) 



The constant depends on the scheme: for a, 



b[/P J = 5/3; for a v , 6{//5 J =0. 

Similarly, for the lattice Ref. |n| prescribes 

lnq 2 M a 2 = r/I, (11) 
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Figure 1. The charm quark's MS mass vs. lat- 
tice spacing squared. Circles (squares) denote 
Method 1 (2). The curve is a linear fit of the 
average (diamonds). Offsets in a 2 are for clarity. 

where /* now comes from the replacement 

q- 2 ^ q- 2 \n(q 2 a 2 ). (12) 

With no constant, this prescription is for the cou- 
pling defined in Ref. (To) , which coincides with ay 
through next-to-leading order. 

When combining the series to form Eq. (|l|), one 
can combine 9b LM and q^ M in the usual way, 



ln( 9 *V) = (/, 



lat 



tVUlat 



t) 



(13) 



provided the constants used to define the I*s are 
compatible. (Otherwise the final q* has problems 
as ma, mj /i — > 0.) Most straightforward, we find, 
is to use ay and to extract TO c h(m c h) directly 
from Eq. (|l|). The resulting q*s are a few GeV 
but somewhat a dependent. 

4. RESULTS 

We have computed the charmonium spectrum 
for (fi, c sw ) = (5.5, 1.69), (5.7, 1.57), (5.9, 1.50), 
and (6.1, 1.40) jl^] . Our (preliminary) results 
for TOch(w c h) with tadpole-improved perturba- 
tion theory are plotted against a 2 in Fig. [IJ The 
error bar is dominated by the unknown two-loop 
correction to Eq. ([!]), estimated to be twice the 
square of the one-loop term. When the analysis is 
repeated without tadpole improvement, but still 
choosing a(q*) as in Sect. 3.2, the data change 



negligibly. The subdominant uncertainty is from 
the Monte Carlo statistics of M 2 qq- 

Extrapolating the average of the two methods 
linearly in a 2 yields 



m ch {m ch 



1.33 ±0.08 GeV. 



(14) 



The error bar now incorporates uncertainty in the 
extrapolation, e.g., extrapolating linearly in a. 
Note that the quoted result neither explicitly cor- 
rects for, nor assigns an error to, quenching, be- 
cause fh(fi) does not run when fj, < m ||. 

A 6% uncertainty for the charm quark's mass 
is twice the 3% quoted for the top quark's mass 
from collider experiments. Alas, without two- 
loop (or nonperturbative) matching for mod =/= 0, 
top standards will be impossible to achieve. 

This work was carried out with A. El-Khadra, 
P. Mackenzie, B. Mertens, and J. Simone. 
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